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Применение К/5 анализа для исследования 
детерминированного хаоса на примере 
возмущенного осциллятора Ван дер Поля 


Исследование базируется на применении В/$ анализа для численного решения уравнения Ван дер Поля ДЛЯ 
изучения фундаментальных свойств сложных экономических систем, которые демонстрируют многомасштабное 
и мультистабильное поведение, а также сосуществование порядка и хаоса в зависимости от изменения 
управляющего параметра. Основной акцент в работ е делается на изучение поведения показателя Херста при 
различных хаотических режимах и горизонтах инвестирования. Изучены возможности В/5 анализа 
для распознавания характера поведения системы. 


Введение 


Целью данного исследования является определение динамики изменения первого 
и второго показателя Херста при переходах между хаотическими и периодическими 
режимами решения уравнения Ван дер Поля, а также изменением горизонта инвестиро- 
вания (укрупнении выборки); определение закономерности такого изменения для исполь- 
зования как дополнительного фактора при идентификации режима функционирования 
исследуемой системы по временному ряду. 

Описание сложной динамики экономических циклов путем выявления регулярных 
и нерегулярных моделей поведения и режимов переключения между различными ди- 
намическими фазами в экономических временных рядах является ключом к улучшению 
экономического прогнозирования. Статистический анализ фондовых рынков и валют- 
ных рынков продемонстрировал перемежающийся характер нелинейных экономических 
временных рядов, которые показывают негауссовское поведение в функции распре- 
деления вероятности изменения цен и степенную зависимость от частоты в спектраль- 
ной плотности. Нелинейные детерминированные модели экономической динамики 
способны моделировать перемежающие временные ряды, которые показывают переход 
от порядка к хаосу, или от слабого хаоса к сильному хаосу, что может объяснить про- 
исхождение и природу перемежаемости, наблюдаемой в экономических системах. В дан- 
ной работе основным методом исследования является метод К/З Херста. 


Осциллятор Ван дер Поля. Основные свойства 


Рассмотрим поведение численного решения возмущенного уравнения Ван дер 
Поля, которое, как показано в работе [1], является модельным примером поведения 
бизнес циклов: 

+ иг -Пх+х=азш(@). 

Данное уравнение — пример двух связанных осцилляторов: внутреннего нелиней- 
ного осциллятора с собственной частотой и внешнего периодического осциллятора с 
движущей частотой ©. В нелинейной системе частота колебаний изменяется от из- 
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менения управляющих параметров, следовательно, в динамической экономической 
модели поведение нелинейных бизнес-циклов зависит от взаимодействия между этими 
двумя частотами под воздействием управляющего параметра. 

При отсутствии внешнего воздействия (а = 0) начало координат (х, = 0, х, =0) 


является единственным решением и образует неустойчивую неподвижную точку; все 
остальные траектории системы стремятся к единственной орбите (предельному циклу), 
окружающей начало координат (Аллайгед, Запег и УогКе, 1996). В присутствии экзоген- 
ного воздействия, уравнение допускает богатое разнообразие периодических и 
апериодических колебаний при изменении параметров д, со и а. Ра|и и ГащегФога (1987) 
привели примеры бифуркационных диаграмм уравнения путем изменения частоты у? 
и амплитуды а, которые показывают режим синхронизации и удвоения периода. Они 
отметили, что симметричность осциллятора Ван дер Поля приводит к сосуще- 
ствованию двух асимметричных аттракторов. При больших амплитудах внешнего 
воздействия много периодических, апериодических и хаотических аттракторов 
сосуществуют вместе. 

Свое исследование мы сосредоточим на влиянии изменения параметра а, при 
постоянных значениях параметра д =1 и м’= 0.45 и начальных условиях (1,1) или (1,1). 


На рис. 1 показаны глобальные бифуркационные диаграммы для уравнения Ван дер 


Поля при симметричных начальных условиях. 
25 25 
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Рисунок 1 — Бифуркационная диаграмма для уравнения Ван дер Поля 
при изменении а от 0,1 до 1,5 для начальных условий а) (1,1) 6) (1—1) 


Наиболее интересным для исследования является промежуток при а = 0,9 доа = 1,1. 
Ему будет посвящена большая часть исследования данного уравнения. 


Применение В/5 анализа 


1 Периодический режим и переход к хаосу 
через слияние аттракторов 


Даный промежуток охватывает значения управляющего параметра а от 0,9831 
до 0,988. Имеется два аттрактора на этом промежутке. 

Рассмотрим поведение показателя Херста вблизи точки бифуркации типа «седло — 
узел». Справа временной ряд периодический, и соответвующий ему показатель Херста 
равняется Н1 = 1,01 (рис. 2), что свидетельствует об абсолютной трендоустойчивости 
ряда, чего и следовало ожидать от периодического временного ряда. Вторым показателем 
Херста будем называть показатель Херста, построенный по точкам, которые следуют 
после «срыва» основного тренда и закачиваются «срывом» второго тренда. 
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Поведение обоих показателей Херста отличается стабильностью на периодичес- 
ком промежутке с незначительным уменьшением для Н/ и более значительным падением 
второго показателя Херста (Н>) при приближении системы к хаотическому режиму 
(а= 0,988) (рис. 2). 
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Рисунок 2 — Показатель Херста для а = 0,984 и а = 0,985 


Рассмотрим изменение показателя Херста в периодическом режиме при укрупне- 
нии горизонта инвестирования в 5, 10, 20 и 60 раз. Если наблюдения рассматриваемой 
системы производятся каждые сутки, то данное разбиение отвечает недели, двум 
неделям, месяцу и кварталу. Подразумевается, что в расчет берутся только рабочие дни: 
в выходные дни экономические показатели не снимаются. Анализ графиков (рис. 3) 
показывает, что укрупнение практически не меняет первый показатель Херста (на 
графике отсутствует) и немного увеличивает второй. Данная тенденция сохраняется 
и для большей степени укрупнения. 
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Рисунок 3 — Изменение второго показателя Херста в зависимости 
от параметра а и степени укрупнения данных 


При достаточно большом укрупнении данных уже невозможно получить инфор- 
мацию о поведении системы и её анализ становится крайне затруднительным. Это 
показано на рис. 4 в): показатель Херста становится значительно меньше 0,5, что говорит 
о сильной антиперсистентности ряда и его нетрендоустойчивости. Напомним, что из- 
начально мы рассматриваем периодическую систему, первый показатель Херста которой 
равен 1. 

Данный анализ показывает необходимость учитывать особенности системы при 
дискретном съеме её показателей. При определенной маленькой частоте съема ин- 
формация о системе теряется, и её анализ К/$ методом становится невозможным. 
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2 Хаотический режим до локальной бифуркации 
типа «седло — узел» (а < аь= 0,98312) 


Проведем аналогичное исследование показателя Херста для данной системы в 
хаотическом режиме до локальной бифуркации типа «седло — узел». 
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Рисунок 4 — Показатель Херста при а = 0,984 и укрупнении данных 
а) в | раз 6) в 20 раз в) в 60 раз 

График В/5 анализа очень похож на график для периодической функции, но 
решение системы не является периодическим. В данном случае показатель Херста 
отображает непериодический цикл, который не имеет частоты; он имеет только 
среднюю частоту, которую и показывает анализ. 

При укрупнении в 60 раз степень хаотичности системы отыгрывает значительную 
роль. Редкая выборка из хаотической детерминированной системы является антипер- 
систентным временным рядом, поведение которого в значительной степени зависит от 
характера и степени хаотичности. 
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Рисунок 5 — Изменение второго показателя Херста в зависимости от параметра а 
и степени укрупнения данных 
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Из рис. 5 четко просматривается тенденция увеличения второго показателя Херста 
при приближении к точке локальной бифуркации типа «седло-узел». Из этого можно 
сделать вывод, что второй показатель Херста увеличивается при приближении к 
порядку и становится меньше при приближении к хаосу. При хаотическом режиме 
влияние разреживания данных (укрупнения горизонта инвестирования) носит менее 
значительный характер на показатель Херста, т.е. применяя метод К/З анализа, еще 
удается узнать о характере поведения системы. 


3 Слияние аттракторов и последующее поведение системы 
(а>ам = 0,98765) 


При удалении от точки слияния аттракторов второй показатель Херста умень- 
шается, переходя в горизонтальную линию уже приа = 0,992. Это означает, что система 
теряет память о периодическом режиме, и остается только средний период системы, 
отображаемый срывом тренда первого показателя Херста. Значение, близкое к 1 говорит 
о детерменированной природе хаотического режима. 

При дальнейшем увеличении параметра а средний период системы увеличивается, 
и временной ряд приобретает более хаотический характер. 

При удалении от точки глобальной бифуркации первые показатели Херста прибли- 
жаются к общему значению, что говорит о стабилизации общих свойств системы и, при- 
нимая во внимание хаотичное поведение системы, отсутствии устойчивых периодических 
ЦИКЛОВ. 

При большом укрупнении К/$ анализ показывает, что система значительно меняет 
свой тип поведения в зависимости от значения параметра а, что не соответствует действи- 
тельности. Поэтому делаем вывод, что данный тип поведения системы является наиболее 
сложным для К/$ анализа, т.к. даже малое изменение в управляющем параметре влияет на 
поведение хаотического ряда (рис. 6). 
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Рисунок 6 — Первые показатели Херста для различных значений а и укрупнений 


4 Переход от хаоса к порядку (а = 1,02 — 1,1) 


Рассмотрим поведение системы при переходе от хаоса обратно к порядку через 
синхронизацию частот. Бифуркационная диаграмма рассматриваемого промежутка 
представлена на рис. 7. 

При изучении второго показателя Херста обратим внимание на его значительное 
падение при переходе от хаоса к порядку и постепенное увеличение при уменьшении 
количества разнородных частот, т.е. при их синхронизации. Это соотносится с 
наблюдением, приведенным ранее, что увеличение второго показателя Херста является 
признаком увеличения упорядоченности системы. 
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Рисунок 7 — Второй показатель Херста при различных значения а 
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Рисунок 8 — Вторые показатели Херста для различных значений а и укрупнений 


При анализе вторых показателей Херста отметим их неизменность при укрупнении 
данных с незначительным уменьшением при 20-кратном укрупнении (рис. 8). Это говорит 
о высокой устойчивости данного режима системы относительно редких выборок. Однако 
при 60-кратном укрупнении анализ системы с помощью В/5 метода, как и для предыду- 
щих периодов, не является показательным: система теряет информацию о своем характере 
поведения. В/З анализ при редких выборках может быть полезен для общей характери- 
стики системы, т.е. её близости к периодическому либо к хаотическому режиму. При боль- 
шом разбросе хвоста графика К/$ анализа большая вероятность, что система близка к пе- 
риодическому режиму и, наоборот, вытянутость графика и показатель Херста, значительно 
больше 0 говорит о хаотическом режиме системы. 


Выводы 


В работе были рассмотрены фундаментальные свойства и проведен В/$ анализ 
уравнения возмущенного осциллятора Ван дер Поля как прототипа модели динамического 
поведения сложных экономических систем. Данная система объединяет в себе свойства 
мультистабильности (сосуществования аттракторов), мультимасштабности и переме- 
жаемости периодических и хаотических режимов. Выбор данной модели для иссле- 
дования был продиктован её относительной математической простотой, богатым разно- 
образием свойств и широким интересом со стороны экономистов. Основное внимание 
данной работы было направлено на изучение свойств показателей Херста для различных 
режимов системы. 

Разреживание данных или, в экономических терминах, увеличение горизонта ин- 
вестирования обычно приводит к падению второго показателя Херста. Однако если система 
приближается к периодическому режиму — разреживание данных изменяет показатель 
Херста синхронизированно (на пропорционально одинаковую величину относительно всех 
укрупнений) и незначительно снижает его. Существует две возможности исчезновения 
второго показателя Херста: 
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Г) он сливается с первым показателем Херста. Данное явление наблюдалось при 
периодическом режиме системы; 

2) его значение становится равным или близким к 0. То есть график в логарифми- 
ческих координатах становится параллельным оси абсцисс. Данное явление наблюда- 
лось при кризисе слияния аттракторов, когда система теряет память о периодическом 
режиме и остается только средний период системы, отображаемый срывом основного 
тренда (в логарифмических координатах) или первого показателя Херста. 

В рамках исследования замечено, что при определенной степени укрупнения 
выборки становится практически невозможно определить характер поведения сис- 
темы. Показатели Херста становятся меньше 0,5 (даже при изначально периоди- 
ческой системе) и единственным фактором, указывающим на периодический или 
хаотический режим, становится разброс точек в логарифмическом масштабе. При 
периодическом режиме он намного больше и напоминает облако с показателями 
Херста, близкими к 0, тогда как при хаотическом цилиндр с показателями Херста 
больше 0, но меньше 0,5. 

В заключение отметим, что характеристики нелинейных динамических свойств 
экономических временных рядов, полученные К/З анализом, могут быть только 
первым шагом, чтобы понять поведение экономических систем. И только при- 
менение комплексного системного подхода может помочь при прогнозировании, 
принятии стратегических решений и планировании долгосрочной политики. 
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А.К. Горайп, К.5. Уеудотоу 

Тре Обе ое В / $ Апа!у$5 Гог Заду оЁ Эаегииис с, СВао$ оп е Ехатре о? 1е Реги фед Уап 
дег Ро! ОзсШа@г 

Те зблау 15 Базе оп Фе и5е ое К / $ апа[уз15 Юг Фе питепса! зоаНоп оё уап 4ег Ро] едиайоп ю зва4у фе 
Капдатепа! ргорегае$ оЁ сотпр[ех есопопис зузбетл5 фа{ ех и пи!Я-зсае ап по!и$аЫе Бебаутог, аз \ме[ аз 
Фе соех1%епсе оРог4ег ап сБаоз, Черепа оп свапее$ оЁ е сопёго| рагатеег. Те ок ассещ оп Фе зву 
оЁ фе НигзЕ ехропепе БеБаутог аё @Шегеп! сБаойс пло4ез ап4 1еуе!$ оЁ шуезитепе. ТЬе роззЮШез оЁ В / 5 
апа[уз15 ортесостйоп оЁ фе Берах1ог оЁе зузбет аге зва4еа. 
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